Transpose-free Multiple Lanczos and Its Application in Padé Approximation
نویسندگان
چکیده
A transpose-free two-sided nonsymmetric Lanczos method is developed for multiple starting vectors on both the left and right. The method is mathematically equivalent to the two-sided methods without look-ahead or deflation steps, but avoids the use of the transpose of the system matrix. The method is applied to the computation of the matrix Padé approximation to a linear dynamical system. The result is a method which can be labeled Transpose-Free Matrix Padé Via Lanczos (TFMPVL). Under certain circumstances, TFMPVL will actually reduce the total number of matrix-vector products needed. It is illustrated with some numerical examples.
منابع مشابه
Transpose-free Matrix Pad e Via Lanczos Method
A transpose-free two-sided nonsymmetric Lanczos method is developed for multiple starting vectors on both the left and right. The method is applied to the computation of the matrix Pad e approximation to a linear dynamical system. The result is a method which can be labeled Transpose-Free Matrix Pad e Via Lanczos. The method is mathematically equivalent to the two-sided methods, but avoids the ...
متن کاملPadé via Lanczos with Multiple Expansion Points
Through considering a large scale single-input, single-output dynamical linear system, we are led to the task of finding a suitable reduced-order model corresponding to that system. Equivalently, we need to successfully approximate the system’s so-called transfer function. To approximate the transfer function the PVL method considers the Lanczos process and its close knit connection to Padé app...
متن کاملError Estimation of the Padé Approximation of Transfer Functions via the Lanczos Process
Krylov subspace based moment matching algorithms, such as PVL (Padé approximation Via the Lanczos process), have emerged as popular tools for efficient analyses of the impulse response in a large linear circuit. In this work, a new derivation of the PVL algorithm is presented from the matrix point of view. This approach simplifies the mathematical theory and derivation of the algorithm. Moreove...
متن کاملA partial Padé-via-Lanczos method for reduced-order modeling
The classical Lanczos process can be used to efficiently generate Padé approximants of the transfer function of a given single-input single-output time-invariant linear dynamical system. Unfortunately, in general, the resulting reduced-order models based on Padé approximation do not preserve the stability, and possibly passivity, of the original linear dynamical system. In this paper, we descri...
متن کاملEigenvalue Computations Based on IDR
The induced dimension reduction (IDR) method, which has been introduced as a transpose-free Krylov space method for solving nonsymmetric linear systems, can also be used to determine approximate eigenvalues of a matrix or operator. The IDR residual polynomials are the products of a residual polynomial constructed by successively appending linear smoothing factors and the residual polynomials of...
متن کامل